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We have studied pseudoscalar correlation functions computed using the overlap operator. Within the accuracy 
of our calculation we find that the quark mass dependence agrees with the prediction of lowest-order Chiral 
Perturbation Theory xPT for quark masses in the range of m ~ ms/2 — 2ms. We present the results of an 
analysis which assumes lowest-order xPT to be valid to extract the low-energy constants E and /p, as well as 
the strange quark mass. Non-perturbative renormalization is implemented via a matching procedure with data 
obtained using Wilson fermions in the Schrodinger functional set-up. We find that the scalar condensate computed 
here agrees with the one obtained previously through a finite-size scaling analysis. 



1. INTRODUCTION 

Testing for spontaneous chiral symmetry 
breaking in QCD by determining whether a scalar 
condensate is developed and computing its size is 
a challenge tailored for lattice field theory. One 
way to extract the scalar condensate is to look at 
the behaviour of the pseudoscalar mass Mp as a 
function of the quark mass m. Lowest-order chi- 
ral perturbation theory (xPT) predicts this be- 
haviour to be linear ||l|,|j , 

{aMpf = BMam , Bm = 4aS//^ (1) 

where S and f^^ are the only constants of the 
lowest-order effective chiral Lagrangian. In the 
full theory they coincide with the quark conden- 
sate and the pseudoscalar decay constant in the 
chiral limit, respectively. 

Within the numerical accuracy of our calcula- 
tion in the quenched approximation (see below) 
it turns out that the linear behaviour of eq. is 
indeed satisfied for a large range of quark masses, 
extending from m ^ ms/2 to to ^ 2tos. However, 
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the linearity of the data is actually surprising, 
since it is observed in a range of quark masses 
where higher-order terms in the quenched chi- 
ral Lagrangian could be relevant. Furthermore, 
at one-loop order one expects logarithmic correc- 
tions to appear [Qj^]. We suspect that the ob- 
served linearity is actually the result of a large 
cancellation between various higher-order effects. 
A detailed discussion into the role of these effects 
will be presented elsewhere In the prelimi- 
nary analysis presented here, we restrict ourselves 
to lowest-order xPT. 

Our assumption that lowest-order xPT gives 
a good description of the mass behaviour re- 
ceives support from the fact that we have a 
completely independent method for determining 
E, which gives consistent results. Indeed, us- 
ing overlap fermions it was shown that the bare 
subtracted scalar condensate could be obtained 
through finite-size scaling methods Q . This anal- 
ysis has now been repeated by other authors |^,^ , 
who have obtained consistent results. Higher- 
order effects might be important for this method 
as well, though here the corrections will be differ- 
ent. Again, we restrict ourselves to lowest-order 



2 



quenched xPT, which is consistent with the nu- 
merical data being fully aware that this as- 
sumption must eventually be checked. 

In order to be useful for phenomenological ap- 
plications, the bare condensate has to be renor- 
malized. Renormalization factors which relate 
matrix elements of local composite operators in 
lattice regularization to a standard continuum re- 
normalization scheme such as MS can be com- 
puted in lattice perturbation theory. However, 
it is well known that perturbation theory in the 
bare coupling go does not converge very well, and 
although the situation can be improved by con- 
sidering so-called "mean-field improvement" |^], 
it is preferable to determine these factors non- 
perturbatively. One method to compute the re- 
normalization factor for the scalar condensate 
computed using the overlap operator has been 
described in a recent publication Q and will be 
briefly reviewed in the following section. 

2. NON-PERTURBATIVE RENORMA- 
LIZATION 

The chiral Ward identities imply that in lattice 
regularizations that preserve chiral symmetry the 
renormalization factor of the scalar density is the 
inverse of that for the quark mass (see e.g. j9|) 

- = (2) 

The non-perturbative renormalization of quark 
masses via an intermediate Schrodinger func- 
tional (SF) scheme has been studied by ALPHA 
pO| , who have computed the relations between 
the running mass in the SF scheme, msF, and 
the renormalization group invariant (RGI) quark 
mass M , as well as the non-perturbative match- 
ing coefficient between ttisf and the bare cur- 
rent quark mass m-w for 0(a) improved Wilson 
fermions. Thus, in order to determine for over- 
lap fermions, it is sufficient to compute the ratio 
of the bare mass in the overlap Lagrangian, rriov, 
and msF- However, the formulation of the SF 
for overlap fermions is not straightforward, ow- 
ing to the fact that the SF boundary conditions 
are incompatible with the Ginsparg- Wilson rela- 
tion. We have therefore employed the following, 
more indirect approach. 



The relation between the bare mass TOov(5o) 
and the RGI quark mass M is given by 



M ^m(5o)™ov(5o)- 
The ratio Af/mov(5o) can be rewritten as 
M M m^ig'o) 



rriovigo) 



= ZM) 



{romov){go)' 



(3) 

(4) 

(5) 



where gQ ^ go in general, and the hadronic ra- 
dius ro ||ll| is used to set the scale. The factor 
has been computed in for a wide range of 
couplings. The ratio M/mov{go) is then obtained 
by determining (roiriov) and {tqiti^) at a reference 
value Xref of some observable, say Xref = (roMp)^. 
Furthermore, the combination Z^{g'Q){rQm^){gQ) 
is a renormalized, dimensionless quantity. We can 
then define the universal factor Um in the contin- 
uum limit as 

UM=lun{Z^ig'o)iromMo)}, (6) 

Since the renormalization of the scalar condensate 
is the inverse of that of the cjuark mass, we define 
the renormalization factor Zg, which relates the 
bare and RGI condensates as 

(ro?Tiov) 



Zsiga) 



Um 



(7) 



At this point it is clear that all reference to the 
bare coupling gQ has disappeared, and the only 
part that retains a dependence on the lattice reg- 
ularization is (roTOov). 

The quantity Um is in fact the RGI quark mass 



corresponding to {r^MpY 



Xref • 



The results 



from ref. are then easily converted into esti- 
mates for Um at several values of x^cf- We stress 
that no additional calculations are required to de- 
termine the universal factor Um, but that it can 
be derived easily from existing results in the lit- 
erature. Furthermore, using the continuum value 
Um ensures that cutoff effects of O(a^) associated 
with the intermediate use of Wilson fermions at 
bare coupling gQ are entirely removed from the 
renormalization condition. As described in de- 
tail in additional cutoff effects in the case of 
non-zero gQ can be quite substantial. 
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As indicated by eq. (g), the renormalization 
of the scalar condensate and the quark mass are 
inversely proportional if chiral symmetry is pre- 
served. If eq. (0) is rewritten as 

t^ML^^, = ^ ) {(^OTOov)(ffo)}:r„f (8) 

it becomes clear that this condition fixes Zg((7o) 
by requiring that the rcnormalizcd quark mass 
reproduces the continuum result for Wilson 
fermions at every non-zero value of go- In other 
words, since the quark mass rriov is required to fix 
Zs, the latter cannot be used to obtain indepen- 
dent predictions for renormalized quark masses 
for overlap fermions. 

This problem can be circumvented after realiz- 
ing that an alternative renormalization condition 
for Zg is provided by the matrix element of the 
pseudoscalar density. The derivation is entirely 
analogous to the previous case. If we use the no- 
tation 



Gp = (0|P|PS) 



(9) 



as a shorthand for the matrix element of the pseu- 
doscalar density P, and the superscripts "RGI", 
"ov" and "w" to distinguish the RGI matrix el- 
ement from the one in the overlap and 0(a) im- 
proved Wilson regularizations, we find 



(^o'GF)(5(,) 



(10) 



The ratio Gf^^jG^^g'^) is given by |l|,|l§ 

= -^p(5o'Mo) —f (1 + Opamq), (11) 



where TOq is the bare subtracted quark mass for 
Wilson fermions. The improvement coefficient 6p 
is known in one- loop perturbation theory [p^ , 
and the combination Zpmgp/M is known from 
ref . . Again one can consider a reference point 
Xref and define a universal factor in the contin- 
uum limit by 

C/p = lim <^ Zp(ffd,Aio)^^ 

x(l + 6pamq)(rgGg)(g^,)j (12) 

J (roAfp)2=x„f, 



Table 1 

Results for Z^. 







eq. ^ 


eq. 






% 


% 


5.85 


1.5736 


1.05(25) 


1.03(15) 




3.0 


1.04(8) 


1.04(8) 




5.0 


0.99(4) 


0.99(6) 


6.00 


1.5736 


0.98(17) 


1.05(14) 




3.0 


1.03(8) 


1.07(8) 




5.0 


1.00(5) 


1.03(6) 



such that the alternative definition for .Zs(5o) is 
obtained from 



2's(go) 



{rlG°^){m) 



(ro J\/p 



(13) 



As in the previous case one can easily convert the 
results of to obtain the continuum factor [/p. 

Our results for (rorriov) and (rpGp) at several 
values of Xrcf were extracted from the same pseu- 
doscalar correlation functions used to determine 
the condensate according to eq. (^. In Table|l]we 
list our estimates for Zg determined at two values 
of the lattice spacing for several reference points. 
As one can see, the results for the two different 
renormalization conditions are entirely consistent 
within errors. In the following we will use the re- 
sults from eq. ( p^ ) for Xrcf — 3.0. We also note 
that the non-perturbative determination of Zg in 
ref. |l^], which is based on a different technique 
p^], agrees with our data at /? = 6.0. 

3. SCALAR CONDENSATE 

We computed pseudoscalar propagators using 
overlap fermions jl^ at two values of the lattice 
spacing. The parameters of the runs are listed 
in Table |2|, a nd further simulation details can be 
found in |l^,|| . Pseudoscalar meson masses were 
extracted by fitting the numerical data obtained 
for the zero momentum pseudoscalar propagator 
to a standard single cosh form. The masses Mp 
are shown as a function of the bare quark mass in 
Figs. |l| and ^ We observe that Mp is consistent 
with being a linear function of the quark mass in 
the whole region. 
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Table 2 

Simulation parameters. 





P 


Vol 




cfgs 


amov 


5.85 


1.6 


10^ X 


24 


50 


0.047-0.188 


6.00 


1.4 


143 X 


24 


48 


0.028-0.133 
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Figure 1. The pseudoscalar mass (aMp)^ as a 
function of the bare quark mass at /? = 5.85. 



Using eq. ([^) an estimate for aE//^ can be ex- 
tracted easily. Combining these results with the 
renormalization factor Zo, at the relevant lattice 
spacing and using ro/a from jl^ to set the scale, 
we obtain for the RGI condensate the following 
estimates: 



S = 



0.0141(6)(ll)GeV3, 
0.0144(6)(ll)GeV3, 



5.85 
6.0. 



(14) 



Here we have used — 128 MeV and ro — 0.5 fm 
to convert into physical units. The first error is 
the statistical uncertainty in the determination of 
aS//^, whereas the second is due to the error in 
the renormalization factor. The results in eq. (|lj) 
can now be compared with the estimate obtained 
from the finite-size scaling (FSS) analysis a± [3 — 
5.85, which is 



S = 0.0138(16)(10)GeV^ /3 = 5.85. 



(15) 



This result is fully consistent with the above 
value. We note, however, that the scale rg en- 
ters as Tq in eq. (|l5|), and that the extraction 
of the FSS condensate in physical units does not 
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Figure 2. The pseudoscalar mass (aAfp)^ as a 
function of the bare quark mass at (3 = 6.0. 



depend at all on the value of f^. Thus, uncertain- 
ties in the lattice scale affect the two methods of 
extracting E in different ways. 

Of course, eq. (|l|) can also be used to compute 
S from data obtained in simulations with 0(a) 
improved Wilson fermions, despite the fact that 
chiral symmetry is explicitly broken in this for- 
mulation. What makes this strategy attractive 
after all is the fact that high statistics data at 
several values of (3 are available, allowing for a 
controlled continuum extrapolation. Combining 
our Wilson data at /9 = 5.85 with the data from 
the ALPHA Collaboration ||l2| , we have extrapo- 
lated tqE//^ to the continuum limit as a function 
of (a/ro)^ (see Fig. ^). The continuum results 
are entirely consistent with our results obtained 
using overlap fermions at the present level of ac- 
curacy. Fig. H also shows a compilation of results 
from various other determinations pO| , pT 15 1. It 
demonstrates that there is not only broad agree- 
ment of the results but also that lattice spacing 
effects with overlap fermions seem to be rather 
small. 



4. STRANGE QUARK MASS 

The behaviour of the dependence of (aMp)^ 
on the quark mass {am) also allows us to calcu- 
late the renormalization group invariant strange 
quark mass Mg- Using again tq to set the scale 
and rriK — 495 MeV to fix the quark mass we find 
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Figure 3. Comparison of different computations 
of the scalar condensate. We also show the re- 
sult from Wilson fermions ||l^ extrapolated to the 
continuum. The double error bars denote the sta- 
tistical and the quadratically combined statistical 
and systematic uncertainties. 

the results shown in Fig. 4 for Ms + M, where M 
is the average of the RGI up and down quark 
masses. The agreement with the results obtained 
using 0(a) improved Wilson fermions in the con- 
tinuum limit is good, although the error bars in 
the overlap case are rather large due to our small 
statistics. 

5. DECAY CONSTANT 

Since the overlap formalism preserves chiral 
symmetry at non-zero lattice spacing, one can de- 
termine the pseudoscalar decay constant /p with- 
out any further renormalization. By contrast, in 
simulations using Wilson fermions the renorma- 
lization factor of the axial current must be 
computed. In our work we have determined /p 
through the PCAC relation, i.e. 



/p 



2^(0|P|PS), 



(16) 



where TOqv is an input parameter in the sim- 
ulation, while Mp and (0|P|PS) arc extracted 
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Figure 4. The renormalization group invariant 
strange quark mass as a function of the lattice 
spacing. 



from pseudoscalar correlation functions. In order 
to compare results from different fermionic dis- 
cretizations we have plotted ro/p versus (rgMp)^ 
in Fig. H One observes rather good agreement 
among the data employing either overlap (this 
work, jisj) or Domain Wall |22| fermions, as well 
as with the Wilson results extrapolated to the 
continuum limit with non-perturbative estimates 
for Za [§3). The fact that the overlap data at 
f3 — 5.85 are already close to the continuum Wil- 
son results suggests that discretization errors for 
this quantity may be small. By contrast, 0(a) 
improved Wilson results at (3 = 5.85 differ sub- 
stantially from their continuum limit. The appar- 
ent smallness of cutoff effects in the overlap case 
is actually the reason for the excellent agreement 
of Zs estimated either through eq. (0) or ([l^). 



6. CONCLUSIONS 

By employing overlap fermions at two values of 
the lattice spacing in the quenched approximation 
we tested whether or not the scalar condensate E 
obtained from a study of the behaviour of (aMp)^ 
as a function of the quark mass is compatible with 
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Figure 5. tq/p versus (roAfp)^ for various dis- 
cretizations. The lines are fits to our data at 
f3 = 5.85 (upper line) and (3 ~ 6.0 (lower line). 



the one previously extracted from a completely 
independent finite-size scaling method. 

We found this test to come out positive at lead- 
ing order in quenched xPT and postpone the dis- 
cussion of higher-order effects to future work |^] . 
We also determined the pseudoscalar decay con- 
stant, as well as the non-perturbatively renormal- 
ized strange quark mass. Within the relatively 
large statistical errors, the agreement with results 
obtained in the continuum limit using 0(a) im- 
proved Wilson fermions the overlap results 
of ref . and also the Domain Wall formulation 
at non-zero lattice spacing is excellent. There 
are strong indications that cutoff effects for over- 
lap fermions are small, but more statistics as well 
as a larger range of lattice spacings is required to 
corroborate these findings. 
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